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Abstract:  This  paper  reptirta  the  olMervationfi  <>f  aky  krightn>‘m  near  the  aiin  at  altitudes 
up  to  25  km,  using  balloons.  From  the  obser\’c«l  brightness  the  particle  density 
is  deduced  as  a  function  of  altitude  and  divided  into  a  component  diffusing  up 
from  below  and  one  settling  from  above.  The  de<luce<l  influx  of  particles  in  the 
size  range  of  0. 1  to  2.8  ntiemns  radius  is  compared  with  rocket  and  other  observations. 

PeaiaMe:  n  sroft  padore  cooduiaioTCfl  gaHHUc  Had.iioaeHHfi  apKocra  nedu 
iid.iH3H  ro.iHua  Mu  RMCorax  Buuie  25  km,  npoBegeHHue  npH  noMotoM  mapoB* 
30!iAOB.  Ha  Ha<>.iioAeHHit  bpkocth  BUBOAHxacb  naoTHOcTb  uacniuu  icaK 
ipyiiKUHn  BMCOTU  If  paage.ifl.tacb  hr  KOMnoHCHTU,  jHiMiyHjHpyioinMe 
CBepxy.  HiJiie.ieHtiMi)  botok  <incTtm  pa;iiiycoM  0- 1-2-8  MHKpoit  cpaBHHBa^CA 
C,iaHilM.MH  Had.1Io;ieHMH  PUKOT  H  CnyTHHKOB. 


1.  Introduction 

Aa  part  of  a  program  to  invefltigatc  the  feasibility  of  coronal  observatiortR 
from  balloons  the  High  Altitude  Observatory  has  flow'ii  an  externally 
occulted  coronagraph  up  to  heights  of  82  000  feet  to  measure  the  angular 
and  the  wavelength  distribution  of  the  daylight  sky.  As  will  be  demonstrated 
in  this  paper,  the  sky  radiance  close  to  the  sun  is  of  interest  not  only  to  the 
instrument  designer,  who  is  concerned  with  the  contrast  between  a  target 
such  as  a  planet  or  an  artificial  satellite  and  the  sky  at  high  altitudes,  but 
also  to  the  astronomer,  who  can  make  inferences  about  the  sedimentation 
rate  of  meteoric  material  into  the  upper  atmosphere.  These  inferences  can 
serve  as  a  valuable  complement  to  the  direct  det<*rminations  of  meteoric 
influx  made  from  rockets  and  satellites. 


2.  Observations 

The  externally  occulted  coronagraph  (fig.  1)  used  in  this  investigation 
was  a  modification  of  the  instniment  <lesignpd  by  Evans  [I]  and  used 
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Schematic  diagram  of  the  Flying  Coronagraph.  The  external  occulting  diak,  150  cm  from  the  obj^tive  lena, 
an  angle  oi  42  minutes  of  arc.  Diffracted  light  from  the  external  disk  aitd  the  objective  aperture  is  removed  by 
nal  disk  and  the  Lyot  apertiue,  respectively.  A  small  periscope  (not  shown)  just  below  the  external  disk  allows 
the  objective  to  examine  the  skylight  at  scattering  angles  O^O  to  57‘'8. 
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by  one  of  us  [2]  for  ground-based  observations  of  the  solar  aureole.  The 
instniment  proved  to  introduce  scattered  tight  of  radiance  no  higher 
than  10**  of  the  mean  radiance  of  the  solar  disk.  This  represented  an 
improvement  of  approximately  a  factor  of  200  over  the  ordinary  coronagraph. 
The  balloon-borne  coronagraph  was  directed  at  the  center  of  the  solar  disk 
to  within  an  accuracy  of  ±  1  minute  of  arc  by  the  guiding  gimbals  developed 
and  used  in  Project  Stratoscope  I  [3,  4,  5,  6].  The  basic  datum  from  the 
instrument  consisted  of  sky  radiances  recorded  in  three  ways:  (1)  phoCo- 
graphioally  at  scattering  angles  O.^O,  20.°7,  3I.°4,  40.°6,  and  57.°8  from  the 
center  of  the  solar  disk  at  A  =>  0.44  n  with  electric  vector  vibrating  perpend¬ 
icular  to  the  scattering  plane,  (2)  photographically  at  scattering  angles 
firom  l.°67  to  2.°8  with  a  low  dispersion  spectrograph  covering  iiO.37^ 
to  0.79^,  (3)  photoelectricaily  at  a  scattering  angle  10.°3  at  0.52  fi.  The 
spectrograph  and  the  photoelectric  photometer  recorded  total  radiance. 
Standa’dization  was  provided  by  a  sunlit,  calibrated,  opal  glass  and  step- 
wedge  combination  which  was  photographed  simultaneously  with  the  spectro¬ 
graph  and  the  O.^’O  to  67.°8  fields.  The  photoelectric  photometer  was 
standardized  frequently  in  flight  with  a  calibrated  opal  illuminated  by  sun¬ 
light.  An  example  of  the  variation  of  sky  radiance  with  height  at  one 
scattering  angle  appears  in  fig.  2.  The  observations  are  compared  with  what 
would  be  expected,  were  the  sky  radiances  measured  in  a  pure  moleeular 
atmosphere  with  three  different  models  of  the  radiation  originating  from 
below.  Both  the  observations  and  the  calculated  values,  taken  from  [7], 
include  the  effects  of  the  variation  of  the  solar  zenith  distance  during  the 
flight.  Attention  is  called  to  the  sharp  approach  of  the  observed  sky  radiance 
toward  the  calculated  values  above  an  altitude  of  65  000  feet. 

3.  Interpretation 

In  order  to  interpret  observations  such  as  those  shown  in  terms  of  the 
size  distribution  of  aerosols  above  the  balloon  at  any  moment  it  is  necessary 
to  solve  the  equation  governing  the  transfer  of  radiation  through  a  turbid 
atmosphere.  To  this  end,  Sekera  [8]  and  Deirmendjian  [9,  10]  have  demon¬ 
strated  the  convenience  of  dividing  the  equation  of  transfer  into  two  simulta¬ 
neous  equations  —  one  describing  the  intensity  due  to  scattering  by  the 
molecules  and  the  other  describing  the  intensity  due  to  scattering  by  the 
aerosols.  The  first,  governing  the  intensity  in  the  pure  molecular  or  Rayieigh 
atmosphere,  has  been  solved  by  Chandrasekhar  [II]  and  tabulated  by 
Coulson,  Dave,  and  Sekera  [7].  By  considering  that,  insofar  as  the  overall 
transfer  of  radiation  is  concerned,  the  Mie  scattering  in  the  upper  atmosphere 


14tf  O.  NEWKIRK,  JR.  AND  J.  A.  EDDY  [Al 

represents  a  small  perturbation  on  the  solution  of  the  Rayleigh  problem,  we 
can  use  their  tabulated  solution.  Subtraction  of  the  radiance  of  the  pure, 
molecular  sky  from  that  observed  then  yields  the  residual  sky,  which  is 

A(tim) 
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Fig.  2.  Photoelectric  observations  of  the  sky  radiance  at  a  scattering  angle  of  10.°3 
at  4  =3  0.82  fi.  The  "hook”  at  82  000  feet  was  produced  by  the  variation  of  solar  zenith 
distance  during  the  floating  period  of  two  ).  ors.  The  observations  are  compared  with 
three  models  of  the  radiation  from  a  pure  molecular  sky. 

attributed  to  scattering  by  the  larger  particles.  By  assuming  that  the  irradia* 
tion  of  Mie  particles  by  the  rest  of  the  atmosphere  and  by  the  ground  is 
negligible  and  that  ft  ft^,  we  can  write  an  approximate  equation 

^  /  Nir,  h)  H  («,  dr.  (1) 

where 

ft~coeZ, 

zenith  distance  of  the  point  in  the  sky. 
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Bm!  radiance  of  the  nky  at  height  h,  in  the  {Kilarization  state,  aa  due  to 
the  scattering  ))y  large  particles, 

-  mean  radiance  of  the  solar  disk  observed  at  the  height  A, 
to^  -  solid  angle  subtended  by  the  sun. 
h  -  height  in  the  atmosphere. 

ii{t,  ^)  =  yiie  scattering  intensity  function  in  the  jiolarization  state, 
fur  a  spherical  drop  of  size  parameter  t.  of  refractive  index  m, 
scattering  through  an  angle  ^  (see  [12]  for  a  complete  definition), 
scattering  angle  counte<l  from  the  direction  of  forward  scattering, 
iX  -2nr/>l  =  8izo  parameter  of  the  droplet, 

X{r,  h)  =  particle  concentration  in  a  column  above  height  h  with  radii  between 
r  and  r  i-dr. 

The  unit  of  N  is  cm-<  per  micron  dr.  The  usual  convention  u  that  /==  1 
refers  to  radiation  with  electric  vector  vibrating  perpendicular  to  the 
scattering  plane  while  j~'2  refers  to  radiation  with  electric  vector,  normal 
to  this  plane.  Calcidatiun  shows  that  the  errors  committed  by  using  the 


Fig.  3.  The  relative  angular  distribution  of  the  light  scattered  by  large  paiiiciea  for 
altitudes  above  50  000  feet.  All  scattering  diagrams  have  been  normalized  at  4  ^  7* 
a  to  relative  intensity  of  6.  The  algebraic  mean  of  the  measurements  also  appears. 
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simplified  equation  in  the  stratosphere  are  of  the  order  of  15  {lercent.  Since 
most  of  the  published  Mie  scattering  coefficients  are  for  m  =  1 .33,  we  have 
made  our  interpretation  in  terms  of  the  scattering  by  .spherical  water  drops. 
All  particles  are  considere<l  to  have  their  radii  between  the  limits  ri  and  rj. 


\ 

♦ 

f' 
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The  size  diHtribution  of  the  scattering  aerosols  is  taken  to  have  the  form 

wh«e  N  =  No  (rfrt)-*,  (2) 


«J-0,  f/r*<l, 
6^6,  flft>  1- 


Junge  [13]  and  his  oollaboratora  [14,  15]  have  found  by  direct  sampling, 
and  Volx  [16],  Newkirk  [2],  Deirmendjian  [9,  10],  and  Volz  and  Bullrich  [17] 
have  inferred  from  optical  measurements  that  such  a  form  adequately 
describee  the  size  distribution  of  atmospheric  aerosols.  In  principle,  the 


Fig.  5.  A  comparison  between  the  particle  size  distribution  at  66  000  feet  determined 
optically  on  10  September  1960  and  on  3  October  1960  and  the  limits  of  particle 
concentration  set  by  Jiuige  and  his  collaborators  by  direct  sampling  techniquea.  The 
optical  determination  is  ambiguous/ for  r  <  0.1  /i. 
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angular  (Hstrihutvjn  of  the  reHulual  aky,  Huch  as  apiiears  in  fig.  3,  should 
determine  the  parameters  r\,  r^,  ra,  and  d  of  the  size  distribution.  In  fact, 
however,  a  certain  ambiguity  is  inherent  in  the  problem.  This  is  demonstrated 
in  fig.  4  in  which  the  angular  variation  of  scattered  light  from  aerosols  with 
different  size  distributions  is  demonstrated.  A  comparison  lietween  the 
observed  angular  distribution  of  scattered  light  and  that  calculated  shows 
that  the  models  displaye<l  could  hardly  be  distinguished  on  the  basis  of  the 
observations.  (Jlearly,  the  observations  allow  nothing  to  b**  said  concerning 
the  size  distribution  below  a  radius  of  approximately  0. 1  /i.  Between  r  0. 1  ft 
and  3.0  ft,  however,  the  size  distribution  is  fairly  well  determined,  and  the 
parameter  is  established  as  3  <()..' 4.  Examination  of  the  wavelength 
distribution  of  the  residual  sky  does  little  to  remove  the  ambiguity  con¬ 
cerning  the  smaller  particles.  In  spite  of  these  difficulties,  the  optical  data 
yield  a  size  distribution  for  particles  in  the  upper  atmosphere  which  is  quite 
consistent  with  that  obtained  by  direct  sampling  (fig. 


HClCHTfkm) 


Fig.  6.  Observed  variation  of  particle  concentration  with  height  for  r  =  0.15  /< 
determined  optically  on  10  September  1960  and  3  Octolx-r  1060  and  by  Junge  and 
collaborators  by  direct  sampling. 
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By  differentiation  of  the  size  diatrihutiona  with  reHi)ect  to  lieight,  «Mie 
can  determine  the  particle  concentration  at  any  |>articular  altitude.  The 
result  of  such  an  analysis  appears  in  fig.  H.  in  which  the  variation  with 
height  of  particles  of  radius  0.15/4  is  shown.  This  particular  ra<iius  for  the 
particle  was  chosen  in  order  to  allow  a  comparison  with  the  direct  sampling 
results.  Quite  clearly,  our  observations  confirm  the  conclusion  of  Juiige  and 
his  collaborators  that  an  aerosol  layer  at  a  height  of  approximately  65  (KM) 
feet  exists.  The  f<»rm  of  such  a  layer  suggests  that  the  majority  of  particles 
in  it  are  of  terrestrial  origin  and  that  the  maximum  is  due  either  to  the 
concentration  of  particles  by  a/lvection  or  by  in  sifii  formation. 


4.  The  net  vertical  flux  of  particles 


The  sharp  drop  of  particle  concentration  with  height  above  the  maximum 
of  the  aerosol  layer  can  be  further  interpreted  under  the  assumption  that 
a  steady  state  prevails  and  only  two  mechanisms  - diffusion  and  sedimenta¬ 
tion— govern  the  flow  of  particles.  Junge  and  his  collalmrators  [14]  write 
the  vertical  flux  of  particles  y  as 


Y  ^  -  n,w, 


dz 


(3) 


where 


»•  =  concentration  of  particles  of  radius  r, 

Wt  -  fall  speed  of  the  particles  of  radius  r, 
n  ■  concentration  of  air  molecules  (numlier  per  cm^), 
r,  =  concentration  of  particles  relative  to  the  concentration  of  air  mole¬ 
cules, 

D  —  diffusion  parameter. 

The  fall  speed  is  that  given  by  the  Stokes-Cunningham  relation  [IS] 
while  the  diffusion  parameter  contains  the  effects  of  both  eddy  and  molecular 
diffusion.  Within  the  layers  of  the  atmosphere  of  concern  in  our  problem, 
eddy  diffusion  is  dominant  [19].  The  solution  of  eq.  (.1)  in  an  atmosphere 
of  scale  height  H  yields 


n»  -  ««o  e 


-IIB- 


I  (r/m 
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w  \ 


1  -  e 


-  /  cr/Bi  111' 


{*) 


where  the  concentration  is  nto  at  z  -  0.  The  effect  of  increasingly  large  values 
of  the  diffusion  constant  D  on  the  variation  of  particle  concentration  with 
height  appears  in  fig.  7,  in  which  the  diffusion  is  considered  to  occur  from 
a  height  66  000  feet.  By  interpreting  the  obwwveil  drop  in  particle  con- 
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centration  with  height  above  the  maximum  of  the  aerosol  layer  aa  due  to 
the  mechanism  of  sedimentation-diifusion,  we  can  determine  that  the 
diiftision  constant  D  from  66  000  to  80  000  feet  lies  in  the  range  2000  om^/sec 
<D<  10  000  cm^/sec  with  the  Iwst  fit  for  D  6000  cm^/sec.  This  value  is 
in  rough  agreement  with  the  Z)  10^  estimated  by  Lettau  [19]  for  temperate 
latitudes  from  models  of  the  representative  zonal  circulation  and  the  lapse 
rate.  The  fact  that  the  slow  decrease  in  concentration  of  the  curve  m  (fig.  7) 
is  not  observed  indicates  that  the  majority  of  the  particles  in  this  height 
range  are  not  of  meteoric  origin. 

However,  the  additional  assumption  that  D  is  constant  above  80  000  feet 
allows  the  calculation  of  the  total  concentration  (in  the  cm^  column)  overhead 
of  particles  with  radius  r  which  have  been  carried  aloft  by  diffusion.  (In 
fact,  the  usual  profile  of  temperature  va  height  above  80  000  feet  would 


Fig.  7.  Variation  of  relative  particle  concentration  with  height  above  the  maximum 
of  the  aerosol  layer  (taken  at  66  000  feet)  to  be  expecte<l  if  sedimentation-diffusion 
equilibrium  applies.  A  constant  inOux  of  meteoric  material  is  assiuned,  and  the  con¬ 
centrations  are  calculated  for  several  values  of  the  diffusion  constant  D.  If  the  only 
source  of  jwrticles  were  the  sctlimentation  of  meteoric  dust  with  influx  rate  y,  the 
concentrations  m  would  be  observed. 
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lead  to  the  expectation  that  D  should  decrease  sUghtly  with  height.)  This 
oalonlation  yields  column  concentrations  for  0.1  /i<r<3/i  which  are  below 
those  require<l  by  the  scattering  observations.  In  other  words,  the  air  above 
80  000  feet  contains  more  particles  in  this  size  range  uhan  can  be  explained 
by  diffusion  from  below.  Tlie  origin  of  the  excess  particles  is  to  be  found  in 
the  sedimentation  of  meteoric  material  downward.  Thus,  for  a  particular 
value  of  D,  the  observations  lead  to  a  value  of  the  meteoric  influx  y.  In 
fig.  8  appears  a  comparison  between  the  meteoric  influx  derived  in  this 
manner  for  0. 1  ^  <  r  <  3.0  /i  and  as  determined  by  various  rocket  experiments. 
The  balloon  measurements  agree  with  the  more  conventionally  determined 
meteor  fluxes  but  are  somewhat  divergent  from  the  “Venus  Fly  Trap"  [20] 


Fig.  8.  ComiMrison  of  meteoric  influx  rates  inferred  in  this  paper  with  influx  rates 
inferred  from  the  rocket  "Venus  Fly  Trap”  [20],  satellites  [21],  visual  and  radar  meteors 
[22],  zodiacal  light  [23]  and  visual  meteors  [24].  The  higher  HAO  curves  refer  to  the 
results  of  10  September  1960;  the  lower  curves,  to  3  October  1900.  For  =  2000cm*/sec 
the  solid  HAO  curves  apply,  while  for  D  =  10  000  cm*/sec  the  dashed  curves  apply. 
The  lower  limit  for  particles  to  remain  in  the  solar  system  against  radiation  pressure 
and  the  upper  limit  for  particles  to  survive  the  passage  through  the  atnoosphera  [25] 

are  indicated. 


154  O.  NEWKIRK,  JK.  AND  J.  A.  EDDY  [lA 

rwiiltH.  The  cauw  of  thw  difference  is  unknown.  Our  analysis  also  seems 
to  aonfiim  the  calculation  of  Whipple  [25],  that  only  particles  with  radii 
greater  than  about  4  ft  are  d«*stroycd  in  passing  through  the  atmosphere. 

This  study  was  carried  out  under  the  auspices  of  the  Office  of  Naval 
Research,  Contract  No.  303(05). 
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